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Abstract 

This is a short review of the statistical mechanical definition of en- 
tropy production for systems composed of a large number of interacting 
components. Emphasis is on open systems driven away from equilibrium 
where the entropy production can be identified with a logarithmic ratio 
of microstate multiplicities of the original macrostate with respect to the 
time-reversed state. A special role is taken by Gibbs measures for the 
stationary spatio-temporal distribution of trajectories. The mean entropy 
production is always non-negative and it is zero only when the system 
is in equilibrium. The fluctuations of the entropy production satisfy a 
symmetry first observed in |7| and then derived in [11 for the phase space 
contraction rate in a class of strongly chaotic dynamical systems. Aspects 
of the general framework are illustrated via a bulk driven diffusive lattice 



1 Introduction 

The production of entropy in spatially extended systems was already discussed, 
both at length and in depth, by the founding fathers of thermodynamics and 
statistical mechanics. It is instructive to divide the discussion in two (very un- 
equal) parts. 

First there is the situation of a perfectly closed (or truly isolated) classical 
mechanical system undergoing a Hamiltonian dynamics and evolving towards 
equilibrium from nonequilibrium initial conditions. That is known as the ap- 
proach or the convergence or the relaxation towards equilibrium. The issue here 
is to understand the usual second law of thermodynamics that associates an ar- 
row of time to the macroscopic behavior leading to equilibrium and to derive 
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from the reversible microscopic laws the irreversible kinetic and hydrodynamic 
equations describing this evolution on the appropriate space-time scale. 
Secondly, there are the many phenomena of dissipation in open driven systems. 
Here one considers a stationary situation of a system or medium in contact with 
several reservoirs that interact with the system. Think of heat baths or particle 
reservoirs at different temperatures or concentrations through which a nonequi- 
librium state is maintained in the system or also of particle systems subject 
to an external driving field. As a result, a stationary current is installed and 
entropy is produced. 

The plan is that I start by recalling some of the main ingredients in the 
qualitative understanding of macroscopic irreversibility for the first scenario, 
return to equilibrium. That has obtained most attention in the past and it 
is well understood. On the other hand, for the nonequilibrium steady state 
scenario the literature is mostly limited to the close to equilibrium situation and 
Section || must remind us of the phenomenology of entropy production. Then, 
I get a chance to say something new when in Section |] I give a Boltzmann- 
like definition of entropy production in a simplified context. Formula (4.2) is 
most important. The general framework is presented in Sections H and |6[ The 
application in the final Section |?] is on the asymmetric exclusion process. 

The paper uses and studies aspects of thermodynamics, statistical mechanics 
and probability theory with a variety of language and tools. One application 
and source of inspiration is however not included, that of the theory of smooth 
dynamical systems. Yet, some of the work that I will present here was started 
in [19 1 (see also |24l |lj|) from trying to understand the Gallavotti-Cohen results 
in [11] on the fluctuations of the phase space contraction in dynamical systems 
with very strong hyperbolicity assumptions and the 'new theoretical ideas in 
nonequilibrium statistical mechanics' a la Ruelle, see p7| . Some apology can be 
read in Section |^. 

The present paper is, except for the Sections ^ and 7.2, a short review and, by 
lack of space, the many details and proofs are omitted. More can be found in 
[DiJ-(23] and I am especially grateful to Frank Redig for the joint work. 



2 Approach to equilibrium 

Here I only remind the reader about the problem of irreversibility in closed 
Hamiltonian systems and about its resolution. There is nothing new here (see 
e.g. |l3|, |l7], |j) except perhaps for some measures of emphasis. 

What is the problem? The microscopic dynamics is time- reversal invariant. 
Why then do we always see individual macroscopic systems taking a particular 
course in time, evolving towards equilibrium? What makes the direction of time 
in macroscopic behavior? 

The resolution of this paradox has many sides and it can easily take you away for 
quite some time. Yet the beginning is easy and it says that there is absolutely 
no problem, it is a matter of counting and you could not expect otherwise: 
equilibrium is what you should expect. Here the basic ingredient is that there 
is a huge difference in scales between the microworld and the macroworld and 
that irreversibility belongs to the macroworld. 
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2.1 Boltzmann entropy 

Take a system of N (large!) components x(i), i = 1, . . . , N, each having n pos- 
sible values. A microstate x — (x(i),i — 1, . .., N) is one of the n N possible 
configurations. The set of all these configurations (or, when additional con- 
straints are given, some subset of this) is the microscopic phase space tt. From 
a microscopic perspective each element of fl is equivalent and we remain indif- 
ferent in the estimates of plausibility: each x has the same probability for being 
realized. From a macroscopic point of view, some distinction can be made. 
For this we make a choice of macroscopic variables which are approximately 
additive (and, in presence of a 'natural' dynamics, are locally conserved). The 
macrovariable partitions 17 in a number of subsets. Say for il = { + 1, — 1}^ and 
macrovariable Xjy(x) = eacn °f these subsets will contain between 

1 and, for N very large, almost 2 N elements. So we all expect that randomly 
picking an element from f2 gives us a configuration that 'typically' has zero Xn 
(with standard error 1/^/N). In other words, while all microstates are equiva- 
lent, some microstates are more typical from a macroscopic point of view. 
To connect the microscopic world with macroscopic behavior Boltzmann intro- 
duced the entropy 

S(X) = \ogW(X) 

where W(X) is the number of microstates x for which Xjy(x) = X. The loga- 
rithm is very convenient and makes S(X) extensive. In the same way (and by 
some abuse of notation), we speak about the entropy S of a microstate x: 

S(x) = S(X N (x)). 

Of course, the phase space for realistic systems is a bit more complicated. We 
should be speaking about a (classical) Hamiltonian dynamics for ./V point parti- 
cles enclosed in a finite box with x corresponding to the positions and momenta 
of all the particles. The phase space f2 is now the collection of all these mi- 
crostates constrained to given values of certain macrovariables Yjv ■ As we have 
a closed system with energy conservation, we can keep in mind that the total 
energy Y/v(x) = E is thus fixed. We are interested in the evolution of a second 
class of macrovariables Xm- Let us think here about the macroscopic variable 
that corresponds to the spatial density profile, say the number of particles to 
the left of our box. Yet, for this more complicated phase space, more or less the 
same counting procedure can be followed. In this same language, microstates 
have a larger entropy when they correspond to a macroscopic value that can 
be microscopically realized in more ways. Now it is important to get a feeling 
of the huge differences that can arise when the system is large. For air, under 
normal circumstances, there are about 10 25 molecules in a box of 1 to 3 . If we 
divide this box in pieces of 1 cm 3 we get 1060000000000000000000000000 possible ar _ 

rangements that would correspond to a homogeneous density profile (about the 
same number of particles in each piece) while only 10 6 possibilities of having all 
of them in just one such piece of the box. When considering these huge numbers 
and reasoning about plausibilities just as we do in everyday life, it appears that 
only a great conspiracy of initial conditions or dynamics can lead to an effective 
decrease in entropy: that is, we expect that S(xt) grows as xt follows the mi- 
croscopic trajectory and equilibrium corresponds to maximal entropy. The law 
of increase of entropy is a statistical law, but one with a 'moral certainty' 
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It must be added that larger entropy does not necessarily mean less ordered or 
more spatially homogeneous. It depends on the type of interaction and/or the 
value of the energy that is fixed. It suffices to think of gravity (where clustering 
of mass is typical) or of molecular interactions at low energy leading to ordered 
structures. 

2.2 Initial conditions 

Note that the counting argument above is symmetric in time, valid as such for 
prediction as well as for for retrodiction. So the real surprise is to see nonequi- 
librium initial conditions. When trying to see the origin of these special low 
entropy initial conditions we soon get in a chain of arguments leading to cosmo- 
logical questions. Here input from the standard model of physical cosmology is 
necessary. The upshot is summarized by Richard Feynman: "it is necessary to 
add to the physical laws the hypothesis that in the past the universe was more 
ordered, in the technical sense, than it is today... to make an understanding of 
the irreversibility" or, much earlier, by Ludwig Boltzmann: "That in nature 
the transition from a probable to an improbable state does not take place as 
often as the converse, can be explained by assuming a very improbable initial 
state of the entire universe surrounding us. This is a reasonable assumption to 
make, since it enables us to explain the facts of experience, and one should not 
expect to be able to deduce it from anything more fundamental." I refer to the 
pleasant Chapter 7 in |2j| for some specific analysis. 

3 Phenomenology of steady state entropy pro- 
duction 

The definition of entropy production for nonequilibrium steady states cannot 
be completely arbitrary. Standard treatments are however largely restricted to 
close to equilibrium phenomenology. 

Entropy production appears in the close to equilibrium thermodynamics 
of steady state irreversible processes as the product of thermodynamic fluxes 
and thermodynamic forces. The forces are produced by gradients in intensive 
variables, the so called affinities that are maintained throughout. One thinks 
of gradients in chemical or electrostatic potential or in temperature. The fluxes 
are currents in the conjugate extensive variables. In linear response, they are 
linear combinations of the affinities. In this way, linear transport coefficients 
are defined that are functions of the intensive variables that locally characterize 
the state of the system. 

As an example take a cylindric solid material through which a stationary heat 
current Jq is maintained by coupling the material at its right and left ends to 
reservoirs at temperature T r and Ti respectively. Assuming a linear relation 
between the affinity Vl/T and the current Jq with a linear response coefficient 
Lq, one arrives at Fourier's law 

Jq = L Q V(^) = -hVT 

where h = Lq/T 2 is the heat conductivity. Here is expressed the empirical 
finding (close to equilibrium) of the proportionality of the heat current with the 
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temperature gradient. The entropy production in the material is then 



S = JQvi = /^ 2 [V(i)] 2 

Of course, the (close to equilibrium) entropy of the material remains constant 
(its macrostate is unchanged) and all produced entropy is carried away by the 
entropy current to be poured into the reservoirs. (The relation between entropy 
current and entropy production is established in the so called entropy balance 
equation.) If we integrate the entropy production over the material, we see that 
the entropy of the reservoirs is changed by 

dS 1 1 

A = aM ¥ r T e ] 

where a is the cross-section area of the cylinder. To maintain the steady state 
we will have to carry away this extra entropy outside the coupled system. I 
refer to [Q for the standard treatment and to ||, |6) for more recent studies. 

4 Multiplicity under constraints 

An elementary mathematical clarification of the connection between Boltzmann 
entropy and thermodynamic equilibrium entropy goes by considering the mul- 
tiplicity of microstates for a particular macroscopic observation, (see also fl^|). 
Take N particles each of which can be in a certain phase space cell i (for ex- 
ample having energy Ei),i = 1, . . . ,n. In total, there are n N microstates. The 
number of such microstates with mi particles in cell i = 1, ... and m n particles 
in cell i = n (mi + . . . + m n = N) is given by the multinomial coefficient 

N\ 

W N (m l7 m 2 , ■ ■ ■ ,m„) 



mi! . . . m„! 



With proportions pt = rrii/N and via Stirling's formula, 

1 ™ 

i=l 

where the right hand side is the Shannon entropy of the probability measure 
(j>i). Imagine now as further constraint that 

n 
i=l 

Then, the multiplicities or the Shannon entropy of (fl.l| ) is maximal when 



Pi = 

Zf3 

with P = (3(E) found from the constraint. This maximal entropy is the equilib- 
rium or Gibbs' entropy 

S(E) =lnZf} + /3E 
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We can now move E to a new value E + dE while supposing that the energy 
levels Ei do not change. Obviously, the maximal entropy does change and it is 
doing so according to 

dS = (3dE 

which is Clausius' formula dS = SQ/T for the change in equilibrium entropy 
from a heat transfer SQ at reservoir temperature T — 1/(3 for fixed volume. 

I suggest a very similar scheme for entropy production. I sketch it here 
without explicit reference to dynamics just to emphasize, in this highly simplified 
setting, the structure of the definition that later, in fuller glory, will follow. (In 
brackets I will announce the analogue.) 

Consider again N variables a(i),i = l...,N with, for simplicity, a(i) = 0, 1, 
i.e., n = 2 in the above (this should be thought of as the pathspace, the space 
of all microscopic space-time trajectories. Thus N is the number of particles 
times the number of times they are observed; a(i) = 1 indicates a right-moving 
particle at that time while cr(i) = indicates left-moving). Let us select the first 
M < N of these variables and let (r)(i),i — 1, . . . , M) be a fixed configuration 
on them. (That is, we select a particular history for a given space-time volume.) 
That of course constrains the number m\ of all variables that are in the first 
microstate (1): 

M M 
N- M + y~] rj(i) >m 1 >^2 r](i) 
i=i i=i 

Let W^\mm) be the number of microstates a G {0, 1}^ so that a(i) = rj(i) 
for i — 1, . . . , M and so that 

N 

— mN 

i=l 

(That will correspond to a given macroscopic value for some current.) Consider 
the involution rj(i) — > fj(i) = 1 — f](i), i = 1, . . . , M (time-reversal will take that 
role later). I suggest now to inspect the logarithmic ratio 

$m = \n— m - (4.2) 

W y M >(r),m) 

(That will be the entropy production.) Evaluating this for large N, while fixing 
all the rest, gives 

M 

S M (r),m) = ]imS%P(ri,m) = V[r?(i) - fj(i)]\n- (4.3) 

n f— ' 1 — m 



(Observe that for A = lnm/(l — m) indeed exp A/(l + exp A) = m so that (4.3) 
is the product of the field A ^ for ra ^ l/2(= the equilibrium value) and the 
variable current or flux in the space-time volume M.) If we now also let M grow 
very large, we will have that 



M M 



^2 V(i) ~ mM,y\j(i) - fj(i)] w (2m - l)M 



i=l 



G 



and hence, 5'm(?? 5 w)/M converges to the relative entropy 

777 777 1 — 777 

S(m) = (2m — 1) In = m In h (1 — m) In 

1 — 777 1 — 777 777 

(That is, the mean entropy production is positive and it equals the relative 
entropy between the original distribution and its time-reversal.) 



5 Gibbs measures with an involution 

The following must be considered as a program, the ideal context that we should 
keep in mind for more concrete realizations. 

r T7d+l 

I write K for a finite set and take f2 = K "-' . This (d + 1)— dimensional 
configuration space announces that f2 plays the role of pathspace; its elements 
g = (<r t (i), (t,i) G Z x Z d ) are space-time trajectories with values a t (i) 6 K at 
the space-time point (t, i). Consider regular space-time cubes Vt.l = {(M) £ 
: \t\ < T,\i\ < L} centered around the origin, where T stipulates the 
temporal and L the spatial extension. Let 0t,l be the involution on O that 
reverses the time: 

{®T,L<r)t(i) = <T-t(i)X(t,i)eVr,L, (5.4) 
= o"t(i) otherwise 

(I could also have included a kinematical time-reversal n as involution on K 
and then write (8j L (r)((j) = n(a-t(i)) for (t,i) £ Vt.l but, for simplicity, I 
will stick here to the choice 7r ^identity.) Given a local function / on Q, it 
is possible to find large enough T ,L so that for all T > T a ,L > L Q ,f(a) = 
f(at(i),(t 1 i) € Vt,l)- There is therefore no ambiguity in writing fQ for the 
new (time- reversed) function and similarly, for a probability measure \x on £7, to 
write /iO for the new probability measure with expectations /i0(/) = [i(f<d). 
Next, consider a translation invariant space-time interaction potential U = (Ua) 
parametrized by the finite subsets A of Each Ua is a function of the 

variables (T t (i), (t, i) £ A, and let me assume, just for convenience, that Ua = 
whenever the diameter of the set A is larger than a finite radius. I define the 
associated entropy production in a finite set A as 



5 A = SP = J2 Pa® - Ua] (5.5) 

AcA 

Clearly, S\ is asymmetric under time-reversal and it is measurable from the 



values of the variables in the set A. Definition ( |5.5| ) is the analogue of (4.2). 
In all realistic realizations of this definition, S\ can be interpreted as a sum 
over products of fields and currents. The reader is probably waiting to see some 
dynamics and it may seem strange to speak about such an entropy production 
but I only deal with the general framework here and one illustration is contained 
in Section |t]. 

Let (i be a translation- invariant probability measure on Q. I define the mean 
entropy production (MEP) in jj, as the expectation 

MEP([/, M ) = lim^ M (S A ) (5.6) 
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Limits are understood in the sense of increasing cubes A = Vt,l- It is important 
to realize that you do not in fact need to take the /i— average in the above but 
for a fixed (large enough) T,S/\A\ will become //—almost surely equal to the 
MEP if fi, is ergodic with respect to spatial translations. 

Let \x be a translation invariant Gibbs measure for the potential (Ua)- It has an 
entropy density s(\i) that equals the entropy density s(/j,Q) of the time-reversed 
Gibbs distribution. Yet, [L and /x0 need not be equal; they can be discriminated 
via their relative entropy density s([i\fi<d) — s(fi<3\/j,) (see fill ). 
Theorem 1 (MEP) Let /i be a translation invariant Gibbs measure for the 
potential (Ua)- Then 

MEP(C7, fi) = s(n\nO) > (5.7) 

with equality if and only if the potentials U and UQ are physically equivalent. 
Remarks: 

• The identification of the MEP with the relative entropy density was an- 
nounced at the very end of Section 

• There can be no spontaneous breaking of time-reversal symmetry: If the 
MEP is zero, then [i and //0 must be Gibbs measures for the same po- 
tential, which means that (Ua) and (Ua&) must be physically equivalent. 
This property of 'no current without heat' is established for infinite inter- 
acting particle systems in [^3], jgjj . 

I now present the local (in space) fluctuation theorem (LFT) for Gibbs mea- 
sures with an involution. Take T and L large and fix the volume V = Vt,l and 
a subset A = Vt,L' with L' << L. For a function / of the variables in V, define 
the expectation 

V aeK v 

with Zy the normalizing partition function. There is no need for translation 
invariance here and the potential (Ua) is allowed to contain time-reversal in- 
variant hard core interactions. I write 



^v\a(°'a) = 



E 



e 



for the partition function in A c = V \ A with boundary condition cta in A. Put 



FX(a A )=ln- 



Zv\a(&a) 



which, for local interactions, depends on the variables in (the interior boundary 
<9A of) A. 

Put Rl = S A -FX. 
Theorem 2 (LFT) 

• For every function G 

1E V [G(-Rl)} = JE v [G(Rl)e- R -] (5.8) 
Ev[G(-5a)] = JE v [G(S A )e-^ +F ^ } (5.9) 
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• For every family (Xa) of complex numbers 



E v [ e -E AcA ^[^e-^]] =]Ev[e -E.cA( 1 - A -)^ e -^ e ^] (5.10) 
Remarks: 

• I take it understood that, for sufficiently local interactions, the potential 
F]( is of the order of the boundary of A. Therefore, asymptotically in 
logarithmic sense, ( pTs| ) leads directly to the symmetry (^) for the distri- 
bution of the entropy production Sa itself, a symmetry first uncovered in 
|7|, [I!]] for the large deviation rate function of the phase space contraction 
in the SRB state of reversible dissipative mixing Anosov diffeomorphisms. 
One should however not take for granted that is uniformly bounded 
by \dA\. In many of the models to which one wishes to apply the above 
scheme, one really needs to prove that it concerns here just a 'boundary 
term.' One origin of problems can be that time is taken continuous and 
an unbounded number of changes can happen in any finite interval. This 
is less of a problem for jump processes where one uses that the Poisson 
process has all exponential moments but when working on non-compact 
phase spaces, the 'boundary term' F]( can possibly have non-existing ex- 
ponential moments and thus really can change the large deviation rate 
function. 



The identity (5.10) generates, by suitable differentiation, equalities be- 



tween correlation functions. A discussion with an application to the On- 



sager reciprocity relations is contained in |19|. I like to compare (5.10) 
with the Ward identities as we know them from quantum field theory; 
the mathematical origin is very similar and non-perturbative, liberating 
us from close to equilibrium assumptions. 



6 Gibbsian hypothesis 

The ambition in the above framework for the study of noncquilibrium driven 
systems is to use the standard Gibbs formalism and in particular its fluctuation 
theory in the 'current '-ensemble. The mathematics will therefore not deviate 
significantly from what is e.g. found in the contributions ||, |2^, [T^| . 
A frequently asked question is where these Gibbs measures come from or how 
they should be connected with a dynamics. Here two short answers. 

6.1 Pathspace measure construction 

The first answer is illustrated below in Section [?] and many more examples can be 
found in |[9]]-|23) and |24|, [l5| ||. One constructs the Gibbs measure explicitly 
as the pathspace measure of a stationary process. Most easy is the case of 
stochastic dynamics. It is not really necessary to realize the pathspace measure 
as a bona fide Gibbs measure as we know it say for lattice spin systems; what is 
needed is to understand how the pathspace measure is governed by a space-time 
action which is approximately local and additive in space-time. Technically, this 
is a matter of setting up the appropriate Girsanov formula and I refer to the 
standard treatments in 0, . Interestingly enough, this method of constructing 
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the process via a Gibbsian space-time distribution has also lead to new existence 
and uniqueness results for classes of diffusion processes, see e.g. [ p6[ . 
Deterministic dynamics are here more of a problem but, for example, we have 
learnt from p[-|ll| how the Gibbsian structure can be exploited on the level of 
the symbolic dynamics for sufficiently strongly chaotic dynamics. 

6.2 Maximum entropy principle 

The first answer has the advantage of being explicit and directly connected to 
model dynamics that one may have in mind. I believe however that the second 
answer is more to the point. The reason to use Gibbs measures here is exactly 
the same as for using Gibbs measures in usual classical or quantum equilibrium 
statistical mechanics. All that changes is the type of ensemble because we 
must work with the currents as macro-observables. To do this on space-time 
is the only price that must be paid but it is a necessary one. Gibbs measures 
appear then as solutions of the maximum entropy principle but that need not 
be restricted to equilibrium conditions. 

7 Asymmetric exclusion process 

I illustrate the previous generalities using one concrete model, that of a bulk 
driven diffusive lattice gas. First about the mean entropy production, to check 
that it coincides with what one should expect. Secondly, about the local entropy 
production fluctuations. 

7.1 MEP for ASEP 

I start with the asymmetric exclusion process (ASEP) on a one-dimensional 
ring {1,2,..., £} (periodic boundary conditions). Each site i of the ring is ei- 
ther occupied by a particle (denoted by r)(i) — 1) or is empty (r)(i) = 0). The 
particle-configuration 7/ is subject to a Markovian particle conserving asymmet- 
ric hopping dynamics with rates 

c(i,i+ 1, r,) = r,(i)(l - V (i + + r,(i + l)(l - r,{i))^- (7.11) 

for changing the configuration from r\ to rf ,l+1 obtained by exchanging the oc- 
cupations at sites i and i + In words, particles hop to nearest neighbor sites 
when there is a vacancy at a rate that depends on the direction. E is now an 
external driving field. We can easily obtain the space-time interaction from the 
standard Girsanov formula for Markov chains. The entropy production is then 
the relative action under time reversal and its expectation in a steady state is 
the MEP. 

The product measure p u with uniform density u € [0, 1] is a stationary (non- 
reversible) measure for this dynamics. If we now consider a trajectory (r]t,t £ 
[-T, T]) of the stationary process in which at a certain time, when the configu- 
ration is 77, a particle hops from site i to i + 1, then the time-reversed trajectory 
shows a particle jumping from i + 1 to i. The contribution of this event to the 
entropy production is therefore 

lnc(i,m,7y)-lnc(M-|-l,?7 M+1 ) = ^[^(i) (1 — 77(i + 1) ) — 77(7+ 1 ) (1 — r/(z) )] (7.12) 
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This jump in the trajectory itself happens with a rate c(i, i + 1, rj) and therefore 
the mean entropy production equals 

MEP(u,E) = J C (0,M)ln e ^°^ Pu (d V ) 

= E u (1 -u) sinh(-) (7.13) 

That is correct: the MEP is the product of the field E with the current 
j (it, E) where the current j(u, E) is the expected net number of particles passing 
through a given bond: 

j(u,E) = J p u (dl])c(i,i+l,r])[r](i)(l-ri(i + l))-r](i + l)(l-T](i))} 

E 

— u(l — u) sinh(— ) 

In quadratic approximation (that is close to equilibrium) j{u, E) w u(l —u)E/2 
and, 

MEP(u,£) « (7.14) 
h c 

which is the dissipated heat through a conductor in an electric field E with 
Ohmic conductivity h c = u(l -u)/2 = p u ([£(0)(l - £(1)) - - C(0))] 2 )/4 

given in terms of the variance of the microscopic current (at E = 0). 

7.2 LFT for ASEP 

I take a (2 + 1)— dimensional set-up. There are (spatial) squares V$ = [-L, L] 2 P\ 
Z 2 and A = \—L', L'] 2 n Z 2 with V < L large, and a continuous time interval 
[-T, T] in which we observe the ASEP with the external field E in the horizontal 



direction. That is the 2-dimensional analogue of 7.1 but now on V$ with periodic 
boundary conditions and hopping rates 

e E/2 e -E/2 
c(hj,v) = —^-r]{i){l - T)(j)) + f?(j)(l - r](i)) 

for a horizontal bond (ij — i + e\) with ei the unit vector in the positive 
horizontal direction, and 

c(i,j,v) = ^[v(i)(^ -V(j)) +VU)( 1 -V{i))) 

for a vertical bond (ij = i±e2)- As stationary measure I take again p u (Bernoulli 
with density u) and (r] s (i),s £ [— T, T],i e Vpj denotes the stationary process. 



I show what becomes of relations ( [p.q ) and (5.9). Let IEy[-] denote the expecta- 

y ■ 



tion with respect to the process in V = [— T, T] x Vq with (pathspace) law Py 



For a function / measurable from A = [—T,T] x Aq, 

K\fQT,L>] = K[f®T,L]=K[f d{P f®J' L) ] (7T5) 
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so that, with 

RX = In 



dPy E \A 



and / = G(R A ), the identity ( 7.15 ) is just (p.8[). R\ is computed from a 
Girsanov formula for the non-Markovian point proces dPy \ A and for this I first 
need to identify the intensities (see |l^, Q| ) . In order to have a Gibbsian structure 
(allowing me to pass to ( |5.9| )) these intensities must be the same in the bulk of 
A as they were in the bulk of V . That is easy to verify from considering the 
conditional expectation of a function g in Ao at time t + S given the past history 
in A : 

TEv[g(r)t+5(i),i € A )|»?.(i),« e [-T,t],ie A ] = (7.16) 
TE^[JE^[g{ Vt+s (i),i G Ao)ji7 s (i), s G [-T,t],ie V ]|»7.(t),« e [-T,t],i G A ] 

The conditional expectation inside is explicit from the Markov process in Vb: 

^v[9{Vt+s(i),i&^o)\Vs(i),sG[-T,t],i&V ] (7.17) 
= g(TH(i),i € A ) + 6L$g(rit(i),i G F ) + 0(<S 2 ) 

with Ly the Markov generator of the process. Since is a sum over all bonds 
with local rates, we see that the process restricted to Ao has the same rates 
except for the boundary of A where a birth and death process is added. As a 
result the Girsanov formula for R\ is indeed 



R A = Sa- F a 



measurable in A with 
S A (rj A )=E f hii)^-Vt(i+ei))-Vt(i+ei)(l-Vt{i))W!(t) (7.18) 

J - T i,i+ ei £A 

where Nl(t) is the number of jumps between i and i + e\ up to time t, and 

= EE/>^ («•> 

+,l(!),n S^^ (() 

where the second sum is over all bonds bi starting at site i G Ao with the other 
end j G Ag and, with bi = (ij), 

Agfa.t) = e Ej '/ 2 ]Ey[l — r} t (j)\r) s (k), k G A , s G [-T,t]] 

and 

^(r?,t) = e E *' 2 E%[THtj)\v.{k),ke A ,s G [-T,*]] 

for £Jj = ±i? if j = i ± ei and £y = if j = i ± ea< The e xpres sion ( 7.18 ) is the 
entropy production, that is field times current in A, and ( 7.19 ) is the boundary 
term. That establishes ([T|). 
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